Magnetic correlations in the Hubbard model on triangular and Kagome lattices 
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In order to study the magnetic properties of frustrated metallic systems, we present Quantum 
Monte Carlo data on the magnetic susceptibility of the Hubbard model on triangular and Kagome 
lattices. We show that the underlying lattice structure is important, and determines the nature and 
the doping dependence of the magnetic fluctuations. In particular, in the doped Kagome case we find 
strong short-range magnetic correlations, which makes the metallic Kagome systems a promising 
field for studies of superconductivity. 
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Frustrated spin systems have received significant at- 
tention because of the possibility of novel magnetic 
ground-states and excitations Q. The triangular spin- 
1/2 Heisenberg model has long-range magnetic order in 
the ground state, while the spin- 1/2 Heisenberg model 
on the Kagome lattice is considered to be disordered. 
The discovery of superconductivity at hK in Na^OoC^ ■ 
J/H2O has generated new interest in frustrated interacting 
systems 2j. Furthermore, recently superconductivity has 
been discovered in /3-pyrochlore osmate KOS2O6 with T c 
of 10A" [3J. This is interesting since the pyrochlore lattice 
structure is a three-dimensional analog of the Kagome 
lattice. 

In cobaltates, cobalt and oxygen ions form a two- 
dimensional triangular network. The hopping matrix ele- 
ment of electrons in the cobalt 3d orbitals is not isotropic, 
and it has been shown that the triangular C0O2 lattice 
consists of four coupled Kagome sublattices 0. Hence, it 
is important to compare the magnetic properties of inter- 
acting systems on triangular and Kagome lattices. The 
electronic properties of the t-J and the Hubbard models 
on the triangular and Kagome lattices have been stud- 
ied using various techniques of many-body physics. The 
triangular t-J model was investigated within the RVB 
framework 0, 0, ano - by using high-temperature ex- 
pansions ||- The triangular Hubbard model was stud- 
ied with the path-integral renormalization-group (RG) 
H, the one- loop RG pj| and the fluctuation-exchange 
(FLEX) [llj approaches. The FLEX method was also 
used for studying the magnetic properties of the Hub- 
bard model on the Kagome lattice [12J. 

In this paper, we compare the nature of the magnetic 
correlations in the Hubbard model on the triangular and 
the Kagome lattices using Quantum Monte Carlo (QMC) 
simulations. We consider the Hubbard model on the 
Kagome lattice to be a simple limiting case to explore 
for new physics due to the underlying orbital structure 
in frustrated interacting systems. The orbital degrees 
of freedom create the possibility for mixing the spin and 
charge channels, and hence of new electronic states. This 
is important because of the general ongoing research ef- 



fort on the transition metal oxides. These are our mo- 
tivations for performing the QMC calculations on the 
Hubbard model on the triangular and the Kagome lat- 
tices. We are particularly interested in the doped cases 
of these models for which there are no exact calculations 
on the magnetic properties. 

In the following, we will see that the nature of the 
magnetic fluctuations on the triangular and the Kagome 
lattices are different. For the triangular lattice, the QMC 
results show that there are strong antiferromagnetic (AF) 
correlations near half-filling at low temperatures, when 
the Coulomb repulsion U is of the order of the bandwidth. 
On the other hand, for weak U, the magnetic correlations 
saturate as T — > 0. In the Kagome lattice, the unit cell 
consists of three atoms and the unit cells form a triangu- 
lar lattice, as seen in Fig. 1(a). Consequently, there are 
three bands of magnetic excitations. Two of these modes 
involve enhanced short-range AF correlations. We find 
that, in the doped Kagome case, the low-frequency short- 
range AF correlations are stronger in comparison to the 
triangular lattice. Hence, we note that it would be use- 
ful to investigate the possibility of superconductivity in 
metallic Kagome systems. 

The Hubbard model is defined by 

H = -t (4r c Jff + hx -) + U ^2 n H n il ~ ^^2 n icr,(l) 

(i,j),<7 i i<7 




FIG. 1: (a) Sketch of the triangular and the Kagome lattices. 
The Kagome lattice is obtained from the triangular lattice by 
removing the sites with the open circles, (b) First Brillouin 
zones of the triangular (black) and the Kagome (red) lattices. 
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where t is the hopping matrix element between the near- 
est neighbor-sites, U is the on-site Coulomb repulsion, 
and fi is the chemical potential. Here, Ci a (c irr ) anni- 
hilates (creates) an electron with spin a at site i, and 
riia- = c^Ciu. In the following, we will take t < and 
consider (n) > 1.0, which is the appropriate case for 
the cobaltates~0 . In obtaining the QMC data presented 
here, the determinantal QMC technique 0] was used. 

For the triangular lattice, the longitudinal magnetic 
susceptibility at frequency u> = is defined by 



x(q) 



JO n 



! (m z (r, +ll r)m 2 (r,)), (2) 



where m z (rj) 



and m z (rj,r) 



e m z (ri)e~ . In this paper, results on \ will be shown 
in units of \t\ . 

The Kagome lattice is a three-band model, since each 
unit cell consists of three sites. Hence, each lattice on 
the Kagome lattice site can be represented by the indices 
(£, d) where I is the unit-cell index and d denotes the 
atomic site in a particular unit cell. The corresponding 
Brillouin zone (BZ) is reduced with respect to the trian- 
gular case as illustrated in Fig. 1(b). For the Kagome 
case, we define the longitudinal magnetic susceptibility 
at lo = as 



where m^(i*i) 



/3 



(m z d (r i+i ,T)m z d ,(ri)), (3) 



(4da) is the an " 



nihilation (creation) operator of an electron with spin 
a at lattice site (i, d) and the summation is performed 
over the unit-cell locations. Diagonalizing the 3x3 ma- 
trix Xdd'(q), we obtain x Q (q) which describes the three 
modes of the magnetic excitations on the Kagome lat- 
tice, which we define as the Kagome magnetic bands. At 
this point, it is useful to note that already for the non- 
interacting (U = 0) case, the triangular and the Kagome 
lattices have different properties. For t < 0, the one- 
electron density of states N(u>) of the triangular lattice 
has a van Hove singularity at (n) = 0.5. For the Kagome 
lattice, there is a (5-function singularity in N(to) at the 
bottom of the band and there are van Hove singularities 
at (n) — 1.16 and 1.51. In addition, N(u> — 0) vanishes at 
(n) = 1.33. These features of N(w) are also reflected in 
the magnetic susceptibilities of the noninteracting case. 

We first present results for the triangular lattice at 
half-filling. Figure 2(a) shows x(q) versus q for U = 4\t\ 
on various size lattices as the temperature is lowered. 
Here, it is seen that x(q) has a broad peak centered 
at the K point of the BZ, and hence the system ex- 
hibits short-range AF correlations. We also note that 
x(q) does not vary significantly with T, in particular for 
0.25|t| < T < 0.17|i|. For comparison, Xo(q) for the non- 
interacting system at T = 0.17|i| is shown by the dotted 
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FIG. 2: Magnetic susceptibility x(q) versus q of the trian- 
gular Hubbard model at zero frequency. Here, q is scanned 
along the path T —* M — > K — > T in the BZ of the triangular 
lattice illustrated in Fig. 1(b). The temperature evolution of 
x(q) at half-filling is shown in (a) for U = 4\t\ and in (b) for 
U = 8|t|. In these figures, the dotted curves represent results 
for the noninteracting case at the lowest temperature used in 
that figure. The evolution of x(q) versus q with the electron 
density (n) is shown in (c) for U = 8|t| and T = . 33 £ | . Here, 
the dotted curve represents the results for the noninteracting 
system at (n) = 1.3 and T = . 33 £ | . 



curve. Figure 2(b) displays x(q) versus q for U = 8\t\ at 
half-filling, where we observe a large Stoner enhancement 
of the AF correlations. In contrast with the U = 4|i| case, 
here, x(q) at the K point grows rapidly with a Curie-like 
T dependence, as T decreases from l\t\ to 0.33|i|. How- 
ever, it is not known whether x(q) saturates at lower T 
for U = 8|i|. These results show that the T dependence 
of x(q) depends strongly on the value of U /\t\ in the tri- 
angular Hubbard model, in agreement with the findings 
of the path- integral RG calculations [fj. 

Figure 2(c) shows the filling dependence of x(q) f° r 
U = 8\t\ while T is kept fixed at 0.33|i|. Here, we ob- 
serve that the AF correlations decay monotonically as 
the electron filling is varied from 1.0 to 1.3. We have 
also performed calculations for x(q) a t higher electron 
fillings. We find that, when (n) is increased to 1.5, the 
peak in x(q) shifts to the M point. Upon further doping 
to (n) = 1.75, we find that, for U = 8\t\ and T = 0.2\t\, 
the Stoner enhancement is about 20%. Hence, for this 
dilute hole concentration, the magnetic correlations are 
weakly affected by the presence of the on-site Coulomb 
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(c) U = 4|t|, (n)=1.15 




FIG. 3: Magnetic susceptibility Xa(<l) versus q for the Hub- 
bard model on the Kagome lattice at zero frequency. Here, 
the three magnetic modes of Xa(q) are shown at each q point, 
as q is scanned along the path F — > M' — * K' — + F in the 
first BZ of the Kagome lattice illustrated in Fig. 1(b). Re- 
sults on Xa(q) at half-filling are shown in (a) for U = 4|t| and 
T = 0.2|f| and in (b) for U = 8\t\ and T = 0.33|t|. Fhese data 
have been obtained on lattices with 6x6 and 4x4 unit cells. 
In (c), results are shown for (n) = 1.15 with U = 4|t| and 
T — 0.14|t|. In these figures, the data points for the 6x6 lat- 
tice are connected by black lines, and the red curves represent 
the results for the noninteracting case. 



repulsion at this temperature. For U — 4|t|, x(q) exhibits 
a slow monotonic decrease with the electron doping away 
from half-filling. 

Next, we discuss the magnetic properties of the Hub- 
bard model on the Kagome lattice. Figure 3(a) shows 
QMC results on x«(q) for U = 4|t| and T = 0.2|t| at 
half- filling. Here, Xa(q) for the three magnetic bands 
are plotted as a function of q, and the red curves repre- 
sent results for the noninteracting case. In this figure, we 
observe that the top band (a = 1) is flat in q space, and 
the second magnetic band (a = 2) is degenerate with the 
first one at the zone center. The third mode is weaker in 
magnitude and exhibits a smooth q dependence. Figure 
3(b) shows Xa(q) versus q for U — 8\t\ and T = 0.33|i| 
at half-filling. The general features are similar to those 
seen in Fig. 3(a), however here the Stoner enhancement 
is larger. In Fig. 3(c), the QMC results are shown for 
U = 4|t| and T = 0.14|t| at (n) = 1.15. These figures 
show that the features of X<*(q) are m correspondence 
with those of the noninteracting case. 
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FIG. 4: (a) xi(ci = 0) versus T for the Kagome lattice with 
U = 4\t\ and 8|t| at (n) = 1.0 and 1.15. (b) Zero-frequency 
component of the magnetic correlations between the nearest- 
neighbor sites C versus T for the Kagome and the triangular 
lattices at (n) = 1.15 for U = 8|£| and 4|i|. In these figures, 
the results on the Kagome lattice were obtained for 6x6 (solid 
points) and 4x4 (open points) unit cells, and the results on 
the triangular lattice are for 12 x 12 (solid points) and 8x8 
(empty points) lattices. In addition, the data points for the 
6x6 Kagome and the 12 x 12 triangular lattices have been 
connected by solid lines. 



In order to gain insight into the origin of the Kagome 
magnetic bands, we note that the Kagome lattice is ob- 
tained from the triangular lattice, as illustrated in Fig. 
1(a), by removing the sites with the empty circles. This is 
equivalent to putting an infinitely repulsive one-electron 
potential at these sites. Bragg scattering from this static 
charge-density-wave field then folds the BZ of the tri- 
angular lattice, and also mixes the different wavevector 
components of the spin fluctuations. It is this process 
which creates the Kagome magnetic bands. 

We next discuss the T dependence of the magnetic 
correlations in the Kagome lattice. Figure 4(a) shows 
Xi(q = 0) versus T for U = 4|*| and 8|i| at (n) = 1.0 
and 1.15. Here, we observe that, at half-filling and for 
U = 4\t\, Xi(q = 0) saturates as T decreases. On 
the other hand, for U = 8|i|, Xi(q = 0) has a strong 
T-dependence. For (n) — 1.15 and U = 4|i|, we ob- 
serve that Xi(q = 0) gets enhanced at low T, and be- 
comes larger than at half-filling. We also see that the 
enhancement of Xi(q = 0) depends on the lattice size. 
However, this type of non-monotonic doping dependence 
was not observed for the triangular Hubbard model. In 
addition, we observe that, for U = 8\t\ and at half- 
filling, xi(q = 0) exhibits a Curie-like T dependence for 
0.25|t| <T< l.0\t\. In the ground state of the Heisen- 
berg model on the Kagome lattice, it is considered that 
a spin gap As w J/20 exists, where J is the magnetic 
exchange Q. For U = 8\t\, we have J « U 2 /U = 0.5|i|, 
which gives A5 « 0.025|t|. We expect that x<*(q) sat- 
urates before T becomes comparable to A<?. However, 
already at T = 0.33|t| we observe enhanced magnetic 
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correlations. 

The FLEX calculations for the Hubbard model on the 
Kagome lattice find that the leading magnetic mode is 
nearly q independent, and the tendency to electronic in- 
stabilities is suppressed We note that while the 
a = 1 mode is nearly q independent, the a = 2 mode 
has ferromagnetic q dependence in the sense that it de- 
creases away from the T point. In addition, the eigen- 
vectors of Xdd'(q) show that, at q « 0, the a — 1 and 2 
modes describe excitations involving the AF polarization 
of the spins within a unit cell |l4| . Away from q = 0, 
these modes have additional structures, however they 
involve AF polarizations over most of the BZ. Hence, 
these two modes contain enhanced short-range AF fluc- 
tuations. Within the context of superconductivity me- 
diated by magnetic fluctuations, an important quantity 
is the zero-frequency component of the magnetic fluctu- 
ations between the two nearest-neighbor sites i and j, 



C= / dr{m z {v u T)m z {v 3 )). (4) 
Jo 

In Fig. 4(b), we compare C versus T for the triangu- 
lar and the Kagome lattices for U — 4|i| and 8\t\ at 
(n) = 1.15. This figure shows that, at these tempera- 
tures, the nearest-neighbor AF correlations are stronger 
for the Kagome lattice, even though the ground-state of 
the spin- 1/2 Heisenberg model has long-range order on 
the triangular lattice and it is disordered in the Kagome 
case. Hence, it would be useful to investigate the possi- 
bility of superconductivity in metallic Kagome systems. 

Finally, we discuss the implications of these data for 
the magnetic correlations observed in the cobaltates. The 
cobaltates have a rich phase diagram with superconduc- 
tivity found for x « 0.35 in Na^OoC^ -y^O |2j and with 
magnetic order and large quasi-particle renormalizations 
observed in Na^OoC^ when x is near 0.75 [l^. For the 
triangular Hubbard model, we find that the magnetic 
correlations are strongest at half-filling, and the correla- 
tion effects are weak in the overdoped regime (n) > 1.5. 
In contrast, for the Kagome lattice we have seen that 
the magnetic correlations can be stronger in the doped 
case when U = 4|i|. However, in both of these models, 
we find that the correlation effects are most prominent 
in the vicinity of half-filling at the temperatures where 
the QMC calculations were performed. We note that it 
would be useful to determine experimentally the q de- 
pendence of the magnetic fluctuations in the supercon- 
ducting Na^CoC^ • yH 2 0. 

In this paper, we have investigated the nature of the 
magnetic correlations in the Hubbard model on the tri- 
angular and Kagome lattices. At the temperatures where 
the QMC calculations were performed, we find, in both of 
these models, that the magnetic correlations grow rapidly 
as T decreases at half- filling for U = 8\t\, while they sat- 



urate when U = 4|t|. In the triangular Hubbard model, 
the AF correlations decay monotonically with the elec- 
tron doping. In the Kagome case, on the other hand, we 
have seen that the magnetic correlations can be stronger 
in the doped case when U = 4\t\. We have also seen that 
in the Kagome case the BZ is reduced and there are three 
modes of magnetic excitations. The two leading modes 
involve short-range AF correlations. In particular, we 
find that the low-frequency short-range AF correlations 
are stronger in the doped Kagome case than in the trian- 
gular case. This makes the interacting metallic systems 
with Kagome type of lattice structures a promising field 
for studies of superconductivity. We conclude that in 
frustrated interacting systems the underlying lattice and 
orbital structures are important in determining the mag- 
netic properties. 
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